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Fig. 1: In GW, the self energy is given by the Hartree term plus a Fock-like term which is 
however in terms of the screened Coulomb interaction W (double wiggled line) instead of the 
bare Coulomb interaction V (single wiggled line). The interacting Green function G is denoted 
by a double straight line. 

1 Introduction 

An alternative to density functional theory HI |2] f3j for calculating materials ab initio is the 
so-called GW approach BE). The name stems from the way the self energy is calculated in 
this approach: It is given by the product of the Green function G and the screened Coulomb 
interaction W, see Fig. [TJ 



This self energy contribution supplements the standard Hartree term (first diagram in Fig. CD) 
Here, r and r' denote two positions in real space, u the frequency of interest; and the imaginary 
unit i in front of G stems from the standard definition of the real time (or real frequency) Green 
function and the rules for evaluating the diagram Fig.[H see e.g. J6). 

In GW, the screened Coulomb interaction W is calculated within the random phase approxi- 
mation (RPA) [7J. That is, the screening is given by the (repeated) interaction with independent 
electron-hole pairs, see Fig.|2] For example, the physical interpretation of the second diagram in 
the second line of Fig. [2]is as follows: two electrons do not interact directly with each other as 
in the first term (bare Coulomb interaction) but the interaction is mediated via a virtual electron- 
hole pair (Green function bubble). Also included are repeated screening processes of this kind 
(third term etc.). Because of the virtual particle-hole pair(s), charge is redistributed dynami- 
cally, and the electrons only see the other electrons through a screened interaction. Since this 
RPA-screening is a very important contribution from the physical point of view and since, at the 
same time, it is difficult to go beyond, one often restricts oneself to this approximation when 
considering screening, as it is done in GW. 

For visualizing what kind of electronic correlations are included in the GW approach, we can 
replace the screened interaction in Fig. Q] by its RPA expansion Fig. |2] This generates the 
Feynman diagrams of Fig. El We see that besides the Hartree and Fock term (first line), addi- 
tional diagrams emerge (second line). By the definition that correlations are what goes beyond 
Hartree-Fock, these diagrams constitute the electronic correlations of the GW approximation. 
They give rise to quasiparticle renormalizations and finite quasiparticle lifetimes as well as to 
renormalizations of the gaps in band insulators or semiconductors. 

It is quite obvious that a restriction of the electronic correlations to only the second line of 
Fig. El is not sufficient if electronic correlations are truly strong such as in transition metal 
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Fig. 2: The screened interaction W (double wiggled line) is calculated from the bare Coulomb 
interaction V (single wiggled line) and corrections which describe screening processes. In 
GW, the screening is given by the random phase approximation, i.e., only bubble diagrams in 
a geometric series (second line) are considered. This geometric series can be generated from 
a single bubble connected to the screened interaction (first line). If we start with W = V on 
the right hand side of the first line, we will generate the second term of the second line, and by 
further iterations obtain the whole series (whole second line). 






Fig. 3: Substituting the screened interaction of Fig.\2\into the self energy diagram Fig.\T\gen- 
erates besides the Hartree term (which is independently considered) the Fock term (first line), 
and some electronic correlations beyond (second line). 



oxides or /-electron systems. The GW approximation cannot describe Hubbard side bands or 
Mott-Hubbard metal-insulator transitions. Its strength is for weakly correlated electron systems 
and, in particular, for semiconductors. For these, extended sp 3 orbitals lead to an important 
contribution of the non-local exchange. This is not well included (underrated) in the local 
exchange-correlation potential V xc of the local density approximation (LDA) and overrated by 
the bare Fock term. The GW exchange is "in between" in magnitude and energy dependence, 
which is also important, as we will see later. If it is, instead, the local correlation-part of 
V xc which needs to be improved upon, as in the aforementioned transition metal oxides or f- 
electron systems, we need to employ dynamical mean field theory (DMFT) [8llH[T0l or similar 
many-body approaches. 

Historically, the GW approach was put forward by Hedin Q as the simplest approximation 
to the so-called Hedin equations. In Section [2l we will derive these Hedin equations from a 
Feynman-diagrammatical point of view. Section [3TT1 shows how GW arises as an approximation 
to the Hedin equations. In Section l3~2l we will briefly present some typical GW results for 
materials, including the aforementioned quasiparticle renormalizations, lifetimes, and band gap 
enhancements. In Section @J the combination of GW and DMFT is summarized. Finally, as a 
prospective outlook, ab initio dynamical vertex approximation (DrA) is introduced in Section 
[5] as a unifying scheme for all that: GW, DMFT and non-local vertex correlations beyond. 
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2 Hedin equations 

In his seminal paper [4], Hedin noted, when deriving the equations bearing his name: "The 
results [i.e., the Hedin equations] are well known to the Green function people". And indeed, 
what is known as the Hedin equations in the bandstructure community are simply the Heisen- 
berg equation of motion for the self energy (also known as Schwinger-Dyson equation) and the 
standard relations between irreducible and reducible vertex, self energy and Green function, po- 
larization operator and screened interaction. While Hedin gave an elementary derivation with 
only second quantization as a prerequisite, we will discuss these from a Feynman diagram- 
matic point of view as the reader/student shall by now be familiar with this technique from the 
previous chapters/lectures of the Summer School. Our point of view gives a complementary 
perspective and sheds some light to the relation to standard many body theory. For Hedin's 
elementary derivation based on functional derivatives see [4] and [0. 

Let us start with the arguably simplest Hedin equation: the well known Dyson equation, Fig. 
HI which connects self energy and Green function: 

G{lV) = G°{11') + G(12)/3S(22 / )G°(2'1') (2) 
= G°(W) + G°(12)/3S(22')G (2 / 1 / ) + G°(13)/3S(33')G (3'2)/3S(22 / )G°(2 / 1 / ) + . . . 

Here, we have introduced a short-hand notation with 1 representing a space-time coordinate 
(i"i,Ti) also subsuming a spin if present; employ Einstein summation convention; G and G° 
denote the interacting and non-interacting (V = 0) Green function^, respectively. Here and in 
the following part of this Section we will consider the Green function in imaginary time with 
Wick rotation t — > —ir, and closely follow the notation of ifTTlR 

In terms of Feynman diagrams, the Dyson equation means that we collect all one-particle irre- 
ducible diagrams, i.e., all diagrams that do not fall apart into two pieces if one Green function 
line is cut, and call this object S. All Feynman diagrams for the interacting Green function are 
then generated simply by connecting the one-particle irreducible building blocks X by Green 
function lines in the Dyson equation (second line of Fig. |4|. This way, no diagram is counted 
twice since all additional diagrams generated by the Dyson equation are one-particle reducible 
and hence not taken into account for a second time. On the other hand all diagrams are gener- 
ated: the irreducible ones are already contained in £ and the reducible ones have by definition 



'Please recall the definition of the Green function with Wick time-ordering operator T for t\ and r{: 

G(ll') = -(Tc(l)c{iy) . (3) 
= -(Tc(l)c(l'y)G( n - t v ) + (rc(l') t c(l))e(r 1 , - n) (4) 

The first term of Eq. © describes the propagation of a particle from 1' to 1 (and the second line the corresponding 
propagation of a hole). Graphically we hence denote G(ll') by a straight line with an arrow from 1' to 1 (The 
reader might note the reverse oder in G(ll'); we usually apply operators from right to left). 

2 Note that in ifTTl summations are defined as /„ dr or in Matsubara frequencies J2 m ; Fourier transfor- 

mations are defined as G(r) = J2 m e ~ lVmT G{v m ). The advantage of this definition is that the equations have then 
the same form in r and v m . Hence, we employ this notation in this Section and in Section 3. In Sections 1 and 
4, the more standard definitions [6] are employed, i.e., dr; (1//3) $Z m . This results in some factors f3 (inverse 
temperatur), which can be ignored if one only wants to understand the equations. 
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Fig. 4: One of the five Hedin equations is the well-known Dyson equation, connecting the 
interacting Green function G (double line), non-interacting Green function G° (single line) and 
self energy E. 

the form of the second line of Fig. |U 

Let us note that the Dyson equation © can be resolved for 



G(ll') = (j[G°) 1 J 



(5) 

(ii') 

with a matrix inversion in the spatial and temporal coordinates. It is invariant under a basis 
transformation from r x to, say, an orbital basis or from time to frequency (note in a momentum- 
frequency basis, the G and £ matrices are diagonal). 

A second equation of the five Hedin equations actually has the same form as the Dyson equation 
but with the Green function and interaction changing their role. It relates the screened Coulomb 
interaction W to the polarization operator P, see Fig. [51 which is a generalization of Fig. [2]to 
arbitrary polarizations. As W we simply define (sum) all Feynman diagrams which connect to 
the left and right side by interactions V. Physically, this means that we consider, besides the 
bare interaction, also all more complicated processes involving additional electrons (screening). 
Similar as for the Dyson equation ©, we collect all Feynman diagrams which do not fall apart 
into a left and a right side by cutting one interaction line V, and call this object P. From P, we 
can generate all diagrams of W connecting left and right side by a geometric series (second and 
third line of Fig. [5]) with a repeated application of P and V (second line of Fig. [5]). As for the 
Dyson equation, we generate all Feynman diagrams (in this case for W) and count none twice 
this way. 

Mathematically, Fig. [5] translates into a second Hedin equation 

Will 1 ; 22') = 1/(11'; 22') + W(ll'; 33')P(3'3; 4'4)V(44'; 22') . (6) 

Note, that in a general basis the two particle objects have four indices: an incoming particle 2' 
and hole 2, and an outgoing particle 1' and hole 1, with possible four different orbital indices. 
In real space two of the indices are identical ri = r^/ and r 2 = r 2 / § 

Next, we turn to the polarization operator P which can be related to a vertex T*. This is the 
standard relation between two particle Green functions (or response functions) and the vertex 



3 This is obvious for the bare Coulomb interaction y(ri, r^; r 2 , r 2 ) = V(ri, r2)<5(ri — r' 1 )<5(r2 — r 2 ) with 
V r (ri,r2) = | ri l r ,| ■ As one can see in Fig. [5] this property is transfered to W for which hence only a 
polarization with two r's needs to be calculated in real space. 
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Fig. 5: A second of the five Hedin equations is an ananlogon of the Dyson equation, Fig. 
E] but for the Coulomb interaction. It relates the screened Coulomb interaction W ( double 
wiggled line) with the bare Coulomb interaction V and the polarization operator P. As the 
Dyson equation can be considered as the defining equation for E, this second Hedin equation 
effectively defines what P is. 
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Fig. 6: A third Hedin equation relates the polarization operator P to two separated Green func- 
tions ("bubble" term) plus vertex (Y*) corrections. This is the standard relation between two 
particle Green functions and fully reducible vertex Y. However since the polarization operator 
cannot include interaction-reducible diagrams Y has to be replaced by Y* (see text). Second 
line: In terms of real space or momentum (but not in an orbital representation) two indices can 
be contracted to a single one (cf. footnote^. 

and represents a third Hedin equation. That is P is given by the simple connection of left and 
right side by two (separated) Green functions plus vertex corrections, see Fig. [6] 

P(W; 22') = /3G(12 / )G(21 / ) + /3G(13)G(3 , l')r*(33 / ; 44')/3G(4'2 / )G'(24) . (7) 

Note, in real space 2 = 2' and 1 = 1' (cf. footnote [3]) so that working with a two index 
object, see second line of Fig.[6l is possible (and was done by Hedin), the inverse temperature j3 
(ks = 1) arises from the rules for Feynman diagrams in imaginary times/frequencies, see ifTTTl. 
Let us keep in mind, that in P or Y* not all Feynman diagrams are included: those diagrams, 
that can be separated into left and right by cutting a single interaction line have to be explicitly 
excluded, see Fig. |7] This is the reason why we put the symbol * to the vertex Y*; indicating 
that some diagrams of the full vertex Y are missing. 

Having introduced the two-particle vertex Y*, the fourth Hedin equation is obtained by relat- 
ing this vertex to another object, the particle-hole irreducible vertex T* h . This relation is the 
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Fig. 7: Left: A Feynman diagram that is not part ofT* (or P) since cutting a single interaction 
line separates the diagram into left and right part. Right: This Feynman diagram is included in 
T* (or P). 
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Fig. 8: A fourth Hedin equation is the Bethe-Salpeter equation between the irreducible r* h and 
particle-hole reducible vertex T*. The particle-hole irreducible vertex T* h collects the Feynman 
diagrams which cannot be separated into left and right part by cutting two Green function lines. 
All diagrams are then generated by the Bethe-Salpeter equation. 



standard Bethe-Salpeter equation. As in the Dyson equation, we define the vertex T* as the 
irreducible vertex T* h plus a geometric series of repetitions of r* h connected by two Green 
functions, see Fig. [8} 

r*(ll'; 22') = r; h (ll'; 22') + r*(ll'; 33 , )/3G(3 / 4)G(4 / 3)r; h (44 , ; 22') (8) 

Here, the particle-hole irreducible vertex collects all Feynman diagrams which cannot be sepa- 
rated by cutting two Green function lines into a left and right (incoming and outgoing) part. The 
Bethe-Salpeter equation then generates all vertex diagrams by connecting the irreducible build- 
ing blocks with two Green function lines, in analogy to the Dyson equation for the one-particle 
irreducible vertex E. 

Since in the polarization operator P (and the corresponding reducible vertex vertex T*) dia- 
grams which connect left and right by only one bare Coulomb interaction line V are however 
excluded, we have to explicitly take out this bare Coulomb interaction V from the vertex V*, 
i.e., we have the standard particle-hole vertex r ph (i.e., all particle-hole irreducible diagrams) 
minus the bare Coulomb interaction V diagram, see Fig. [9] 

r; h (ll'; 22') = r ph (ll'; 22') - 1/(11'; 22') (9) 

The bare V and any combinations of T* h and V are then generated in the screening equation 
©, Fig. El 
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Fig. 9: We have to exclude the bare Coulomb interaction V from the particle-hole vertex T p h 
since such contributions are already considered in Fig. \5\ 
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Fig. 10: By hands of selected Feynman diagrams we illustrate that differentiation ofT, w.r.t. G 
yields the particle-hole irreducible vertex. 



In this fourth equation, Hedin directly expresses r* h as the derivative of the self-energy w.r.t. 
the Green function Q (respectively V [4]). This is a standard quantum field theoretical relation 



rp h (ll';22') 



5G{2'2) ' 



(10) 



which in terms of Feynman diagrams follows from the observation that differentiation w.r.t. G 
means removing one Green function line, see Fig. [10] and Ref. ifTTll . If we, as Hedin, consider 
the self energy without Hartree term £ H artree(ll') = v ( ll 'i 22')G(2'2) we obtain the vertex T* h 
instead of r ph : 



r; h (ll';22 



(ID 



5G(2'2) 

Note, the derivative of the Hartree term w.r.t. G(2'2) yields V(ll'; 22'), i.e., precisely the term 
not included in T* h , see first diagram of Fig. [Kjl 

The fifth Hedin equation is the Heisenberg equation of motion for the self energy, which 
follows from the derivative of the Green function © w.r.t. n, i.e., the time-part of the coordinate 
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1. Let us start with the Heisenberg equation for the Heisenberg operator c(l) (H = 1): 

9c(ri,ri) 



071 

and a general Hamiltonian of the form 



[c^t^H] (12) 



H = Fo(ll')c(l)t C (l') + -7(11'; 22')c(l)t c (l')c(2)t c (2 / ). (13) 



From Eq. (fT2l) . we obtain the Heisenberg equation of motion for the Green function^) 

ftn(ii'\ 

= 5(1 - 1') + tfo^'M^l') - 7(13'; 22') (rc(3')c(2) t c(2 , )c(l') t ) . (14) 



9ri 

The last term on the right hand side of Eq. (U~4"l) is by definition the self energy times the Green 
function. Hence this combination equals the interaction 7 times a two-particle Green function. 
The two particle Green function in turn is, in analogy to Fig. [6l given by the bubble term (two 
Green function lines; there are actually two terms of this: crossed and non-crossed) and vertex 
correction with the full (reducible) vertex T. Besides the Hartree term, the self energy is given 
by (see Fig. flTT): 

E(ll') = -7(13'; 22')/3G'(4 , 2)G'(2 , 4)G'(3'3)^(31 , ; 44') - 7(12'; 21')G(2'2) (15) 

Note Eq. (fTBl) can be formualted in an alternative way (see Appendix [A] for a detailed calcula- 
tion): Instead of expressing this correlation part of the self energy by the bare interaction and 
the full vertex we can take out the bare interaction line 7 and any particle-hole repetitions of 7 
from the vertex, i.e., take T* instead of T. Consequently we need to replace 7 by W to generate 
the same set of all Feynman diagrams: 

S(ll') = -W{13'; 22 , )/3G(4'2)G(2 / 4)G(3 / 3)r* (31'; 44') - 17(12'; 2l')G(2'2) (16) 

see second line of Fig. [Tlfl 

The five equations Eq. ©, ©, ©, ©, (© correspond to Hedin's equations (A13), (A20), 
(A24), (A22), (A23), respectively 01 [or to equations (44), (46), (38), (45), and (43), respec- 
tively, in [5]). This set of equations is exact; it is equivalent to the text book quantum field 



4 Note, the first term here is generated by the time derivative of the Wick time ordering operator, the second 
term stems from [c(l), H (22')c(2)^ c(2')] and the third one from 

[c(l),iy(33';22')c(3)t c (3') C (2)tc(2')] 

employing [A,BC] = B[A,C] + [A,B]C and the Fermi algebra {c(l'), c(l)t} = c(l')c(l) t + c(l)t c (l') = 
5(l-l'),{c(l) )C (l')} = {c(l)t )C (l')t} = 0. 

5 Note, Hedin defines as a "vertex" A [5 1 a combination of T* and two Green function lines. Because he works 
in real space (or momentum) coordinates only at a common point 2 = 2' needs to be considered (in the same way 
as in the second line of Fig. |6j. Hedin also adds a "1" in form of two ^-functions: 

A(ll'; 2) = r*(ll'; 33')^G(3'2)G(32) + 5(1' - 2)6(2 - 1) . (17) 
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Fig. 11: The fifth Hedin equation is the Heisenberg equation of motion, which connects the 
one-particle Green function with the two particle Green function or as shown in the figure (and 
as employed) the self energy with the vertex. 

theory j6l [EH relations between r,r ir , £ and G; but it contains additional equations since W 
and P are introduced. The advantage is that, this way, one can develop much more directly 
approximations where the screened Coulomb interaction plays a pronounced role such as in the 
GW approach. 

3 GW approximation 

3.1 From Hedin equations to GW 

The simplest approximation is to neglect the vertex corrections completely, i.e., to set set r* h = 
0§ Then the Bethe-Salpeter equation © yields 

T* = 0. (18) 

The polarization in Eq. CD) simplifies to the bubble 

P GW (11'; 22') = /3G(12')G(21') . (19) 

The screened interaction in Eq. © is calculated with this simple polarization 

W{ll'] 22') = 1/(11'; 22') + W(W; 33')P GW (3'3; 4'4)V(44'; 22') . (20) 

The self energy in the Heisenberg equation of motion (PT6l) simplifies to Fig. [2l i.e., 

E GW (11') = -W{12'; 21')G(2'2) (21) 

(plus Hartree term). 

From this, the Green function is obtained via the Dyson equation (0: 

G(ll') = G°(ll') + G(12)/3S GW (22 / )G°(2 / 1 / ) (22) 
These five (self-consistent) equations constitute the GW approximation. 
6 Note this violates the Pauli principle, see last paragraph of Appendix lAl 
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3.2 GW band gaps and quasiparticles 

While the five GW equations above are meant to be solved self-consistently, most calculations 
hitherto started from a LDA bandstructure calculation^ and calculated from the LDA polariza- 
tion (or dielectric constant) a screened interaction W which in turn was used to determine the 
self energy with the Green function Go from the LDA: £ = iGoWo. 

Such calculations are already pretty reliable for semiconductor band gaps, which are underes- 
timated in the LDA. Due to the energy(frequency)-dependence of S bands at different energies 
are under the influence of differently strong screened exchange contributions. In semiconduc- 
tors, it turns out that the conduction band is shifted upwards in an approximately rigid way. 
The valence band is much less affected so that the GW band gap increases in comparison to 
the LDA gap. This effect can be mimicked by a so-called scissors operator, defined as cutting 
the density functional theory (DFT) bandstructure between valence and conduction band and 
moving the conduction band upwards. Cutting LDA bandstructures by a pair of scissors and 
rearranging them yields the GW bandstructure within an error of 0.1 eV for Si and 0.2 eV for 
GaAs KT71. 

More recently, self-consistent GW calculations became possible. Many of these calculations 
employ an approximation of Schilfgaarde and Kotani lfT8l [T9l where instead of the frequency 
dependent GW self energy S nn / (u, q) a frequency-independent Hermitian operator 

S nn / = Re[E nn /(e g , q) + S„/ n (e ? , q)]/2 (23) 

is constructed in the basis n, n' employed in the GWfLDA algorithm. This self energy operator 
has the advantage that (as in LDA) we can remain in a one-particle description and employ the 
Kohn-Sham equations with Hermitian operator S nn / to recalculate electron densities and Bloch 
eigenfunctions. 

The band gaps of this self-consistent approach are slightly larger than experiment, see open 
triangles of Fig. [[2] This can be improved upon and band gaps can be calculated very reliably if 
additional to GW (some) vertex corrections are taken into account. In Fig. [[21 the inclusion of 
electron-hole ladder diagrams (visualized on the right hand side of Fig. [12]) results in the filled 
triangles with band gaps being within a few percent of the experimental ones. As in other areas 
of many-body theory, doing the self-consistency without including vertex corrections does not 
seem to be an improvement w.r.t. the non-self-consistent GW since self-consistency and vertex 
corrections compensate each other in part. Full GW calculations beyond the Schilfgaarde and 
Kotani one-particle-ization (|23T) have only been started and applied to simple systems such as 
molecules lQ"4l[P3l and simple elements lfT5l[T6ll . 

Besides this big success to overcome a severe LDA/DFT shortcoming for semiconductor gaps, 
GW or GoWq calculations also show a quasiparticle renormalization of the bandwidth. For 
alkali metals, electronic correlation are expected to be weak. Nonetheless experiments observe 
e.g. in Na a band narrowing (of the occupied bands) of 0.6 eV [|22l compared to the nearly free 
electron theory. While GW If23l yields such a band narrowing, it is quantitatively with 0.3 eV 

7 An alternative, in particular for / electron systems, is to use LDA+U as a starting point ifTZl . 
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Fig. 12: Theoretical vs. experimental band gap ofDFT, self-consistent GW (scGW) without 
and with (some) vertex corrections (electron-hole ladder diagrams) (reproduced from /f2i]/). 

only half as large as in experiment ll2~3l . Of the more strongly correlated transition metals, Ni 
is best studied: here, the occupied d bandwidth is 1.2 eV smaller in experiment than LDA and 
there is a famous satellite peak at -6eV in the spectrum [24J. While GW 11251 yields a band- 
narrowing of 1 eV which is surprisingly good (see Fig. [T3l) . the satellite is missing. In fact, it 
can be identified as a (lower) Hubbard band whose description requires the inclusion of strong 
local correlation. This is possible by DMFT; and indeed the satellite is found in LDA+DMFT 
11261 and GW+DMFT calculations see next section. 

Besides the mentioned band-narrowing which is associated with a reduced quasiparticle weight 
or effective mass enhancement (related to the real part of the self energy), there is also the 
imaginary part of the GW self energy, which corresponds to a scattering rate. For Ag the GW 
scattering rate is reported to be in close agreement with the experimental one obtained from 
two-photon photoemission [281 . 

4 GW+DMFT 

Since GW yields bandstructures similar to LDA (with the improvements for semiconductors 
discussed in the previous section) substituting the LDA part in LDA+DMFT by GW is very 
appealing from a theoretical point of view: Both approaches GW and DMFT are formulated 
in the same many-body framework, which does not only has the advantage of a more elegant 
combination, but also overcomes two fundamental problems of LDA+DMFT: (i) The screened 
Coulomb interaction employed for d-d or /-/ interactions in DMFT can be straight forwardly 
calculated via W; one does not need an additional constrained LDA approximation ll2"9ll30ll3Tl 
to this end; (ii) the double counting problem, i.e., to subtract the LDA/DFT contribution of the 
local d-d or /-/ interaction which is included a second time in DMFT, can be addressed in a 
rigorous manner since for GW+DMFT we actually know which Feynman diagram is counted 
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Fig. 13: Experimental, LDA and GW bandstructure ofNi (reproduced from /l25l/). 




Fig. 14: The DMFT self energy is calculated from the local contribution of all (one-particle 
irreducible) Feynman diagrams. 

twice. 

Biermann et al. 11271 proposed GW+DMFT, which they discuss from a functional integral point 
of view: a GW functional and a local impurity functional are added; the derivatives yield the 
mixed GW+DMFT equations. From a Feynman diagrammatic point of view, this corresponds 
to adding the GW self energy, Fig. Q] and the DMFT self energy which is just given by the 
local contribution of all (one-particle irreducible) Feynman diagrams, see Fig.QjJ From these, 
the local screened exchange GW and the doubly counted Hartree term need obviously to be 
explicitly subtracted for not counting any diagram twice. 

This results in the algorithm Fig. [T5l Here, we leave the short-hand notation of Section |2] and 
13. II with 1, 2 since GW is diagonal in to and k and DMFT is diagonal in u and site indices. Let 
us briefly discuss the GW+DMFT algorithm step-by-step; for more details see 11321 : 

• In most GW calculations, the starting point is a conventional LDA calculation (or another 
suitably chosen generalized Kohn-Sham calculation), yielding an electron density p(r), 
bandstructure e LDA (k) and also an LDA Green function Gk(w) (the latter is calculated as 
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Do LDA calculation, yielding G k (uS) = [wl+//l-e LDA (fc)] _1 . 




Calculate GW polarization P gw {uj) = -2i J ^G(u + uf)G(uf) . 




If DMFT polarization _p DMFT is known (after the 1st iteration), include it 






P GWDMFJ (k,oj) = P GW (k,u;) - J- / d\P GW (k,oj) + P DMFT H . 

VbzJ 




With this polarization, calculate the screened interaction: 

W(k-uj) = V ee (k)[l - VccWPfau)]- 1 . 




Calculate ^^=J ^G q {r = Q-)W(k - q, 0) and 




Calculate T, GW (r, r'; u) = i J %£-G(r, r'; u + u')W(r, r'; to') . 




Calculate the DMFT self-energy £ DMFT and polarization p DMFT as follows: 






From the local Green function G and old self-energy £ DMFT calculate 






(g°)- 1 (Lo) = G^ 1 (uj) + T, DMVr (uj) £ DMFT =0 in 1st iteration. 






Extract the local screening contributions from W: 

U(u) = [W-\u;) - P DMFT (a;)]- 1 . 






With U and Q , solve impurity problem with effective action 






A=J2 «#L)~VZn+£ / dr^r(rW(r) [/ ira (r-V)^*(r')^(T') ) 

ua I'm Imcru' 

resulting in G and susceptibility X- 






From G and X , calculate £ DMFT (u;) = (^W) - G _1 (u;), 

PDMFT (w) = 17-1(0;) _ p _ [/xt/]" 1 (cj). 




Combine this to the total GW self-energy: 

V GmDMFT (k,u) = V GW (k, u)-Jd\ S GW (fc^) + S Hartree (fc^ S» artree + S DMFT (^). 




From this and G°, calculate G^, ew (w) _1 = G^a;)" 1 - E fc (w). 


Iterate with G k = G k ew until convergence, i.e. \\G k - G£ ew <e. 



Fig. 15: Flow diagram of the GW+DMFT algorithm (reproduced from / l32l/ ). 

in the first line of Fig. \\5\ where bold symbols denote an (orbital) matrix representation). 

• From this Green function, the independent particle polarization operator P GW is calcu- 
lated convoluting two Green functions (2nd line of flow diagram Fig. [15]). Note there is a 
factor of 2 for the spin. 

• From the polarization operator in turn, the local polarization has to be subtracted since 
this can (and has to) be calculated more precisely within DMFT, which includes more 
than the RPA bubble diagram (after the first DMFT iteration). 

• Next, the screened interaction W is calculated from the bare Coulomb interaction V ee and 
the overall polarization operator p GW + DMFT , 

• Now, we are in the position to calculate the GW self energy. The first term is the Hartree 
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diagram, which can be calculated straight forwardly in imaginary time r, yielding £ Hartree 
and the corresponding local contribution £j c artree , which we need to subtract later to avoid 
a double counting as it is also contained in the DMFT. 

• The second diagram is the exchange from Fig. [TJ which has the form G times W for the 
GW self energy. 

• This GW self energy has to be supplemented by the local DMFT self energy, which 
together with the DMFT polarization operator is calculated in the following four steps: 

1. The non-interacting Green function Q° which defines a corresponding Anderson 
impurity model is calculated. 

2. The local (screened) Coulomb interaction U (u) has to be determined without the 
local screening contribution, since the local screening will be again included in the 
DMFT. That is we have to "unscreen" W for these contributions. 

3. The Anderson impurity model defined by Q° and U(u) has to be solved for its 
interacting Green function G(u) and two-particle charge susceptibility \- This is 
numerically certainly the most demanding step. 

4. From this G(u) and Q°, we obtain a new DMFT self energy £(a>) and from the 
charge susceptibility a new DMFT polarization operator. 

• All three terms of the self-energy have now to be added; and the local screened exchange 
and Hartree contribution need to be subtracted to avoid a double counting. 

• From this GW+DMFT self energy we can finally recalculate the Green function and 
iterate until convergence. 

The flow diagram already shows that the GW+DMFT approach is much more involved than 
LDA+DMFT. However, it has the advantage that the double counting problem is solved and 
also the Coulomb interaction is calculated ab initio in a well defined and controlled way. Hence, 
no ad hoc formulas or parameters need to be introduced or adjusted. 

For defining a well defined interface between GW and DMFT a particular problem is that GW 
is naturally formulated in real or k space and is presently implemented, e.g., in the LMTO J3 
or PAW basis EUl . However, on the DMFT side we do need to identify the interacting local d- 
or /-orbitals on the sites of the transition metal or rare earth/lanthanoid sites, respectively. The 
switching between these two representations is non-trivial. It can be done by a downfolding 
il33ll34l or a projection onto Wannier orbitals, e.g., using maximally localized Wannier orbitals 
ll35l [371 or a simpler projection onto the d (or /) part of the wave function within the atomic 
spheres ll36l [38l . However, not only the one-particle wave functions and dispersion relation 
need to be projected onto the interacting subspace but also the interaction itself. To approach 
the latter, a constrained random phase approximation (cRPA) method has been proposed 11391 171 
and improved by disentangling the d(or /)-bands [|40l . The latter improvement now actually 
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Fig. 16: Spectrum (k-integrated) of Ni [left: LDA+DMFT (solid lines), spinpolarized LDA 
(LSDA ) ( dotted lines); right: GW+DMFT]. The two lines represent the minority- and majority- 
spin spectrum respectively. At roughly -6eV, a satellite peak is clearly visible in the majority- 
spin spectrum (reproduced from [26] and K27\l . respectively). 

allows us to do cRPA in practice. For the calculation of the two-particle polarization operators 
and interactions, Aryasetiawan et al. iflTTl even proposed to use a third basis: the optimal product 
basis. 

On the DMFT side, the biggest open challenges are to actually perform the DMFT calcula- 
tions with a frequency dependent Coulomb interaction U (u) and to calculate the DMFT charge 
susceptibility or polarization operator. 

As the fully self-consistent GW+DMFT scheme is a formidable task, Biermann et al. G71 
employed a simplified implementation for their GW+DMFT calculation of Ni, which is ac- 
tually the only GW+DMFT calculation hitherto: For the DMFT impurity problem, only the 
local Coulomb interaction between d orbitals was included and its frequency dependence was 
neglected W(w) ~ W(0). Moreover, only one iteration step has been done, calculating the 
inter- site part of the self energy by GW with the LDA Green function as an input and the intra- 
site part of the self energy by DMFT (with the usual DMFT self-consistency loop). The GW 
polarization operator P GW was calculated from the LDA instead of the GW Green function. 
This is, actually, common practice even for conventional GW calculations which are often of 
the G(Wo form (see Section EOl) . 

Fig. US (right panel) shows the GW^+DMFT k-integrated spectral function of Ni which is similar 
to LDA+DMFT results (left panel). Both approaches yield a satellite peak at ta —6 eV. 

5 All of that: ab intito DrA 

From the Hedin equations, it seems to be much more natural to connect (i) the GW physics 
of screened exchange and (ii) strong, local correlations on the two-particle level than on the 
one particle level as done in Gty+DMFT. In the Hedin equations, the natural starting point is 
the two-particle (particle-hole) irreducible vertex. A generalization of DMFT to n-particle cor- 
relation functions is the dynamical vertex approximation (DrA) [42] which approximates the 
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all fully irreducible 
local diagrams 



i i 



Fig. 17: In ab initio DTA we take as the fully irreducible vertex the bare Coulomb interaction 
plus all local vertex corrections (only one such diagram is exemplarily shown) / !??!/ . 

n-particle fully irreducible^ vertex r ir by the corresponding local contribution of all Feynman 
diagrams. For n — 1 the one-particle irreducible vertex is the self energy so that DrA yields 
the DMFT. For n = 2, we obtain non-local correlations on all length scales and can calculate, 
e.g., the critical exponents of the Hubbard model [|43l . 

Recently, some of us have proposed to use this DrA ab initio for materials calculation [44). The 
fully irreducible vertex r ir is then given by the bare Coulomb interaction, which possibly is non- 
local, and all higher order local Feynman diagrams, see Fig. [FT] From T ir the full (reducible) 
vertex is calculated via the parqet equations [fTTI . The calculation of the local part of r ir only 
requires us to calculate the two-particle Green functions of a single-site Anderson impurity 
model, which is well doable even for realistic multi-orbital models. For the parquet equations 
on the other hand, there has been some recent progress B31 . 

As a simplified version of ab initio DTA one can restrict oneself to a subset of the three channels 
of the parquet equations, as was done in Il42ll43l . In this case one has to solve the Bethe-Salpeter 
equation with the particle-hole irreducible vertex r p h instead of the parquet equations with the 
fully irreducible vertex r ir . That is, our approximation to the Hedin equations is to take the 
local r* h (all Feynman diagrams given by the local Green function and interaction) in the Hedin 
equation ©. In practice, one solves an Anderson impurity model numerically to calculate T* h . 
Full and simplified version of ab initio DTA contain the diagrams (and physics) of GW, DMFT 
as well as non-local correlations which are responsible for (para-)magnons, (quantum) critical- 
ly and "all that". 

Support of the Austrian Science Fund (FWF) through 1597 (Austrian part of FOR 1346 with the 
Deutsche Forschungsgemeinschaft as lead agency) is gratefully acknowledged. 



Appendices 

A Additional steps: equation of motion 

In this appendix a detailed explanation is given how to derive the Hedin equation of motion for 
S, i.e., equation (fT6l) . from the standard equation of motion ([TBI) . In a first step T is expressed 

8 Fully irreducible means, cutting any two Green function lines does not separate the diagram into two parts. It 
is even more restrictive (less diagrams) than the particle-hole irreducible vertex (whose diagrams can be reducible 
e.g. in the particle-particle channel). 
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Fig. 18: Relation between V and V* in terms of Feynman diagrams. 

in terms of T*. In order to keep the notation simple, the arguments of all functions are omitted 
and the functions are considered as operators. 

The starting point of the calculations are the Bethe-Salpeter equations for T and T*: 



r = r ph + rf3GGr ph =► r = r ph (i - pggt^y 1 
r* = r; h + v ph (3GGT* ==> r; h = r*(i + pggt*)- 1 . 

Using equation (PTOl) . i.e., r ph = r* h + V, one gets: 

T*(l + f3GGT*y 1 + V ' 



(24) 



1 - + /3GGT*)- 1 + V 



(25) 



Multiplying this equation with 1 = (1+/3GGT*) 1 (1+PGGT*) and using the standard relations 
for operators, A^B^ 1 = (BAy 1 leads to: 

T= fr(l + /3GCT*)- 1 + V^M/fl-( /3G'G + PGGT*PGG) v) (l+pGGT*\ (26) 

1 p 

where 1 = V~ l V was inserted and the definition for the polarization operator, equation ©. 
Now one can use the second Hedin equation ©, which can be rewritten as W = V(l — PVy 1 . 
Inserting this relation into the equation for T, one arrives at the following result: 



r = r*(i + pGGT*y v~ l w 1 + pggt* )+w[i + pggt 



(27) 



(28) 



Another formulation of the second Hedin equation gives V 1 = W 1 + P = W 1 + (1 + 
(3GGT*)f3GG. Replacing V^ 1 by this expression gives: 

r = r* + T*PGGW(l + PGGT*) + W{1 + PGGT*) 
= T* + T*PGGW + T*PGGWPGGT* + WPGGT* + W . 

This equation shows how the full T is related to the T*. Diagrammaticaily this relation is shown 
in Fig. M 

In the next step T as given in equation (T28l) is inserted into equation (fT5l) . yielding 

S = -GVPGGY - GV 
= -GVj PGG + PGGT*PGG) W - GV({ §GG + pGGT*pGG) W + l)pGGT* - GV . 



p 



p 



(29) 
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From the second Hedin equation it follows that VPW = W — V. Inserting this relation into 
Equation (|29l) yields: 



(30) 



E = -G(W -V)- G(W - V)f3GGT* - GV/3GGT* - GV = 
= -GW - GW(3GGT*, 

which is exactly equation (fT6l) . 

Let us also, at this point, mention that T should satisfy an important relation: 

r(ll';22') = -r(12';21') (31) 

This relation is known as crossing symmetry (see e.g. ifTTI . equation 7.5) and is simply a 
consequence of the Pauli-principle: Exchanging two identical fermions leads to a — sign in the 
wave function. The screened interaction W, however, does not fulfill this crossing symmetry. 
Therefore, setting Y* = 0, as it is done in the GW-approximation leads to T = W (see Fig. IT8l) 
which violates this crossing symmetry, i.e. it violates the Pauli-principle. 
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